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Introduction 

The 3-dimensional incompressible Euler equation ( "3D Euler" ) is as follows: 

(| + u ■ V«) u = -V xP (x G R 3 , t > 0) 

V x ■ u = (x e R 3 , t > 0) 

u(x, 0) = (x G R 3 ) 

with m° a given, smooth, divergence-free, rapidly decreasing vector field on R 3 . Here, 
u(x, t) and p(x, t) are the unknown velocity and pressure for an ideal, incompressible fluid 
flow at zero viscosity An outstanding open problem is to determine whether a 3D Euler 
solution can develop a singularity at a finite time T. A classic result of Beale-Kato-Majda 
[1] asserts that, if a singularity forms at time T, then the vorticity u(x, t) = V x x u(x, t) 
grows so rapidly that 

T 

/ sup | u(x, t) | dt = oo. 

J X 



In [2], Constantin-Fefferman-Majda showed that, if the velocity remain bounded up to 
the time T of singularity formation, then the vorticity direction u)(x,t)/ \u(x,t)\ cannot 
remain uniformly Lipschitz continuous up to time T. 

One scenario for possible formation of a singularity in a 3D Euler solution is a con- 
stricting vortex tube. Recall that a vortex line in a fluid is an arc on an integral curve 
of the vorticity u(x,t) for fixed t, and a vortex tube is a tubular neighborhood in R 3 
arising as a union of vortex lines. In numerical simulations of 3D Euler solutions, one 
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routinely sees that vortex tubes grow longer and thinner, while bending and twisting. 
If the thickness of a piece of a vortex tube becomes zero in finite time, then one has a 
singular solution of 3D Euler. It is not known whether this can happen. 

Our purpose here is to adapt our work [3], [4] on two-dimensional flows to three 
dimensions, for application to 3D Euler. We introduce below the notion of a "regular 
tube" . Under the mild assumption that 

T 

\ sup \u(x, t) \dt < oo, 

J x 


we show that a regular tube cannot reach zero thickness at time T. In particular, for 3D 
Euler solutions, a vortex tube cannot reach zero thickness in finite time, unless it bends 
and twists so violently that no part of it forms a regular tube. This significantly sharpens 
the conclusion of [2] for possible singularities of 3D Euler solutions arising from vortex 
tubes. On the other hand, [2] applies to arbitrary singularities of 3D Euler solutions, 
while our results apply to "regular tubes" . 

Although we are mainly interested in 3D Euler solutions, our result is stated for 
arbitrary incompressible flows in 3 dimensions. The proof is simple and elementary. The 
main novelty for readers familiar with [3], [4] is that we can adapt the ideas of [3] to three 
dimensions, even though there is no scalar that plays the role of the stream function on 
M 2 . 

1 Regular Tubes 

Let Q = h x I 2 x J 3 C M 3 be a closed rectangular box (with Ij a bounded interval), and 
let T > be given. 

A regular tube is an open set Q t C Q parametrized by time t G [0, T), having the 
form 

Q t = {(xi,X2,x 3 ) G Q : 9(x 1 ,x 2 ,x 3 ,t) < 0} (1) 

with 

9eC\Qx[0,T)), (2) 
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and satisfying the following properties: 



^x U x 2 6 I 7^ for (x u x 2 , x 3 , t) eQ x [0, T), 9(x 1 ,x 2 , x 3 , t) = ; 



(3) 



^(^3) := {{xi,x 2 ) e h x 7 2 : (^1,^2,^3) G ^} is non-empty, 



(4) 



for all x 3 El 3 ,te [0,T); 



closure (fi t (x3)) C interior (7i x 7 2 ) 



(5) 



for all x 3 el 3 ,te [0,T). 

For example, a thin tubular neighborhood of a curve T forms a regular tube, provided 
the tangent vector V stays transverse to the (xi,x 2 ) plane. 

Let u(x, t) = (uk(x, t))i<k<3 be a C 1 velocity field defined on Q x [0, T). We say that 
the regular tube Q t moves with the velocity field u, if we have 



It is well-known that a vortex tube arising from a 3D Euler solution moves with the fluid 
velocity. 



Theorem: Let Q t C Q(t G [0,T)) be a regular tube that moves with a C 1 , divergence 
free velocity field u(x, i). 



( 



d_ 
di 



+ u . V x = whenever (x, t) E Q x [0, T), 9(x, t) = . 



(6) 



2 Statement of the Main Result 



if 



r 




(7) 



then 



liminf Vol(Q t ) > 0. 



(8) 
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3 Calculus Formulas for Regular Tubes 

Let £l t be a regular tube, as in (1)- • -(5). Recall that 

MtM = {(x u x 2 ) e Ji x J 2 : 0(x 1: x 2 ,x 3 ,t) < 0}. (9) 

Define also 

StM = {(x u x 2 ) G interior (h x I 2 ) : 9(x 1 ,x 2 ,x 3 ,t) = 0} for x 3 G I 3 , t G [0,T). (fO) 

Also, for intervals I C I 3 , and for t G [0,T), define 

fi t (7) = {(xi,x 2 ,x 3 ) eQ : x 3 e I and 0(zi, x 2 , x 3 , t) < 0}, and (fl) 
S t (I) = {(x 1 ,x 2 ,x 3 ) G Q : x 3 E I and (xi,x 2 ) G ^(xa)}. (12) 

Let i/ denote the outward-pointing unit normal to S t (I 3 ), and let v = (i>!,i> 2 ,0), where 
(z>i,i> 2 ) is the outward-pointing unit normal to S t (x 3 ). Thus, v and v are continuous 
vector- valued functions, defined on S = {(xi,x 2 ,x 3 ,t) G Qx[0,T) : (xi,x 2 ,X3) G S t (I 3 )}. 
Define also scalar-valued functions a, a on S by requiring that 

(J^ + av = ^ + ai> • = for x G S t (I 3 ). (13) 

Again, cr and a are well-defined and continuous on S, thanks to (3). Let F be any 
continuous function on Q. We will establish the following elementary formulas. 



d 
It 



J F(xi,x 2 ,x 3 ) d(Area) = J F(x 1 , x 2 , x 3 ) a(x 1 , x 2 , x 3 , t) d (length) 



(xi,X2)eSt(x3) 

(14) 



J Frf(Area) = J | J -^d (length) 



F 

i 

St(I) x 3 €l {St(x 3 ) 



> dx 3 . (15) 



To check (14) and (15), we may use a partition of unity to reduce to the case in which F 
is supported in a small neighborhood U . Also, we may restrict attention to a small time 
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interval J. In a small enough U x J, we may assume that S t (I 3 ) is given by the graph of 
a C 1 function ip, thanks to (3). Thus, without loss of generality, we may suppose that, 
in U x J, we have 



(xi,x 2 ,x 3 ) G Qt(h) if and only if x\ < ip(x 2 , x 3 , t), and 
(x 1 ,x 2 ,x 3 ) G S t (I 3 ) if and only if xi = ijj(x 2 ,x 3 ,t). 

From (16), (17) we have also 

(xi,x 2 ) G £lt( x 3) if and only if x\ < ip(x 2 ,x 3 ,t), and 
(x 1 ,x 2 ) G S t (x 3 ) if and only if x 1 = ip(x 2 ,x 3 ,t) . 



(16) 
(17) 



(18) 
(19) 



In view of (16) • • • (19), we have 



d , dip 



dx 2 dx 3 







d 



1 + — V +h^-0 and 



(20) 



1+ (A^ on5n([/x7). 



(21) 



From (10), (17), we have 9(ip(x 2: x 3 , t), x 2 , x 3 , t) — on S D (U x J). Differentiating in 
t, we obtain 

[^+f^(x 2 ,X3,0,0,oVv a; N ) ^ = on 5 n ({/ x J). 



Subtracting this from (13), we find that (#2, £3, t), 0, 0) — cri/ is orthogonal to V x 6, 
hence also to v. This yields the formula 
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Similarly, subtracting the two equations (13), we learn that {av — au) is orthogonal to 
V x 9, hence also to v. Therefore, 



a = a {y • v) on S . 



(23) 



From (20) • • • (23), we obtain 



a 



' dip s 
~dt 



1 + 



dip 



on «S n (U x J) . 



(24) 



Now we can read off (14) and (15). In fact, (18) gives 



d_ 

dt 



J Fd (Area) 



d_ 
It 



Fdxi dx 2 



j F(ip,x 2 ,x 3 ) — dx 2 , 



and (19), (24) yield 



St 0*3 ) 



J Fad (length) = J Fa ■ ^ 1 + (J^^J dx 2 = J F(if>, x 2 , x 3 ) ^ dx 2 , 



proving (14). Similarly, (17) gives 



J Fd (Area) = J F(ip,x 2 ,x 3 ) ■ + (^"^) dx2 dxs 

St(I) 



while (20) and (21) imply 



v ■ v 



n i+ (£*)/n 



so that (19) yields 
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X3&1 



\.S t (x 3 ) 



V ■ V 



— d (length) 



> dx? 



-/ / 



F 



dx2 



> dx? 



' d_ 

dx-. 



' _d_ 

dx? 



tp dx2 dx? 



This completes the proof of (15). Note that (23) allows us to rewrite (15) in the form 



J Frf(Area) = J < J 



Stir) 



{St(x 3 ) 



F-d (length) 



> dx?, 



(25) 



4 Proof of the Theorem 



We retain the notation of the previous sections. We will define a time-dependent interval 



J t = [A(t),B(t)]ch 



(26) 



and establish an obvious formula for the time derivative of Vol VL t {Jt)- We assume that 
the endpoints A{t), B(i) are C 1 functions of t. We have 



VolQ t (J t ) = J Area Q t (x 3 )dx 3 , so that 



Xz&Jt 



d 
It 



Vo\n t (J t ) = B'(t) Area tt t (B(t)) - A'(t) Areatt t (A(t)) + f — Areafi t (x 3 ) dx 3 . 



Applying (14) with F = 1, we find that 
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d 
~dt 



Vol Q t (J t ) = B'(t) Area fl t (B(t)) -A'(t) Area Q t (A(t))+ J j J ad(length) 



x 3 €Jt [St(x 3 ) 



In view of (25) (with F = a on S), this is equivalent to 

j t Vol n t (J t ) = B'(t) Area Q t {B(t)) - A'(t) Area fi t (A(i)) + ^ <jd(Area). (27) 

St(Jt) 

Now we bring in the hypothesis that fit moves with a divergence-free C 1 velocity field u. 
From (6) and (13), we see that {av — u) - V x 9 = on S t (J t ). Thus (cru — u) is orthogonal 
to v, so that a — u ■ v on S t (J t ), and (27) may be rewritten as 

d 

VolUi(Jt) = Area S2 t - A'{t) Area il t {A{t)) -f 

St(Jt) 



j t Vol Slt(Jt) = B'{t) Area Vt t {B{t)) - A'(t) Area fi t (A(t)) + ^ u • i/d (Area) . (28) 



On the other hand, since u = (ui,u 2 , u 3 ) is divergence-free, the divergence theorem yields 

0= J (V x • u) d(Vol) = y u ■ z/d(Area) + / -u 3 <i(Area) — / -u 3 (i (Area) . 
n*(Jt) St(Jt) n t (B(t)) n t (A(t)) 



Hence, (28) may be rewritten in the form 
d 

voni t (j t ) = 

n t (B(t)) n t (A(t)) 



j t Vo\n t (J t )= J [B'(t)-u 3 (x,t)}d(Aie8,) - J [A'(t)-u 3 (x,t)]d(Area). (29) 



This is our final formula for the time derivative of Vol Qt(Jt)- It is intuitively clear. 

We now pick the time-dependent interval Jt = [A(t), B(t)] C I3. Let I 3 = [a, b], and let 
to £ (0, T) be a time to be picked below. We define 



1 

B(t) = b— / max \u(x, r)| dr , (30) 
t 
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and 



i 

A(t) — a+ / max \u(x, r)| dr . (31) 
J xeQ 
t 

T 

We are assuming that u(x, r) is continuous on Q x [0, T), and that J max \u(x, r)\dr < oo. 

o 

It follows that A(t),B(t) are C 1 functions on [0, T), and that 



a < A(t) < B(t) < b for t E [t , T) , (32) 

provided we pick to close enough to T. We pick to so that (32) holds. Thus, flt(Jt) C Q 
for t G [to,T). Immediately from (30), (31), we obtain 

B'(t) = -A'(t) = max \u(x,t)\ > max \u 3 (x,t)\ (33) 

(recall m = (tti, tt2, u^)). 

From (29) and (33) we see at once that 

^VolQ t (J t ) > forte [t ,T). (34) 
at 

On the other hand, (4) and (32) show that VolQ to (J to ) > 0. 
Consequently, 



limmf VolQ t > limmf Vo\Q t ( J t ) > VolQ to (J to ) > 0. 



The proof of our theorem is complete. 
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